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Neutron Compton Scattering measurements in a variety of materials have shown a relative deficit 
in the total signal from hydrogen compared to deuterium and heavier ions. We show here that a 
breakdown in the Born-Oppenheimer approximation in the final states of the scattering process 
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5h ' Deep Inelastic Neutron Scattering experiments on a variety of systems have found the integrated intensity of 
^ the signal from hydrogen to be smaller than expected, relative to heavy ions such as oxygen, or even deuterium, 
[1] given the well known nuclear cross-sections for these materials. This deficit has also been observed in electron 
Compton scattering. [2] This has been attributed to previously unknown physics at the very short time scales of 
the measurements [3], quantum entanglement between the struck proton and the other protons of the material [4], 
inadequacies in the treatment of the data [5,6], and a failure of the usual van Hove scattering formula to describe 
the data [7] . We show here that the effect can be explained within the usual formalism by a breakdown of the Born 
Oppenheimer approximation in the final state of the scattering process, which contains a very rapidly moving proton 
^ , with sufficient energy to mix the electronic states of the system. These corrections are present for any scattering 
Q • experiment, but only become significant for high momentum transfer. 

O ' We begin with the usual expression for the scattering cross-section when the spin of the neutron is not observed, 
in the large q limit in which only incoherent scattering remains. 

^ : Siq,Lo) = < 0|e-«~^Mi?« > I'Siu; - (1) 

iy~j , where the ground state is taken to be the Born Oppenheimer ground state, 

m ■ . ^ ^ 

O ■ < r, R\0 >= $o.o(i?)ao(r, R) (2) 

^ : 

' r are the coordinates of the electrons and R the coordinates of the ions. The \Efi > are energy eigenstates of the 
complete system of electrons and ions. aQ{f, R) is the solution of the electronic problem with the protons and other 
ions at the positions R, and ^ofiiR) is the ground state for the heavy particles in the potential energy surface defined 
I ' by \ao >. i?i is the component of R denoting the position of the struck particle. There are in addition, a complete 
^ ,: set of excited electronic states, for the ions at the same positions, denoted by ai(r, R). It is evident from Eq. Ithat 

O: r 

U- S{q,uj)du^l (3) 



k>( ' for any q. The total intensity of the incoherent scattering is proportional to this integral times the total scattering 
■ cross-section for the scattering particle. Our point of view is that if this appears not to be the case, then some 
' intensity is not being accounted for in the experimental integration. 

It is usually assumed at this point, that the energy eigenstates of the full Hamiltonian can be adequately represented 
by the ground state for the electronic system and the excited states for the ions in the Born Oppenheimer potential 
for that ground state. The energies Efi are the sum of the electronic energy, E^{R) and the energy of the heavy ions 
in the effective potential defined by this energy surface. However, we wish to consider the limit that the transferred 
momentum is large, so that the recoiling proton has sufficient energy to mix the electronic states. In short, we must 
go beyond the Born Oppenheimer approximation in expanding the state e^''-^^ |0 > in energy eigenstates, and consider 
the full Hamiltonian for the system, which includes the kinetic energy of the struck proton. We will treat this kinetic 
energy as a perturbation, to the extent that it acts on the electronic wavefunctions, and first calculate the energy 
eigenstates to lowest order. We take the full Hamiltonian to be 
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H = H^{r,R) + ^,^ (4) 
The BO energy eigenstates ,|$n,j > |an >satisfy 

Ho (f, E)a„ (r , ^) = < (f, ^) (5) 

and 

(Si^ + = Enj^nj{R) (6) 

In addition to the kinetic energies of the heavy particles, obtained from (Sj2]^|$ >)|q!„ > the leading term in the 
perturbation when acting on a state of the form |<& > |a„ > is 

EiPi|$ > .Pi\an > /M (7) 
The action of the momentum of the heavy particle on the electronic states is, from Eq. 5 

< am\Pi\an >= in < am\-^\an > /(K - ^m) (8) 

6 Hi 

where the matrix elements are taken over the electronic coordinates, and the positions of the heavy particles are the 

same for all the states. This perturbation mixes the original Born-Oppenheimer energy eigenstates. For all the heavy 
particles except the one struck by the neutron, the momentum operator acting on the wavefunction of the particle 
gives a small overall contribution, which is the basis for the usual Born-Oppenheimer approximation. However for 
the struck particle, particle 1 in Eq. 1, this operator gives a large term, of order Hq, and cannot be neglected. 
We have then that the energy eigenstates for the full Hamiltonian are given to first order by 

\En,i >= \^n,i > \a„ > +^mj < ^mj\ < "mlTT; ^ , E^l"" > •TtI*".^ > l*'".^ > I"™ > (9) 

We will neglect for the moment the shift of the energies for these states, as it is not essential for this work. The 
important point, is that because of the presence of the fast proton, the excited states have a significant overlap with 
the ground state, and as a consequence, there is some amplitude for the scattered neutron to transfer some of its 
energy to the electronic system. When Eq. 9 is used in Eq. 1, we see that the only component of the excited states 
wave function that is significant is that which couples it to the ground state. The sum in Eq. 1 extends over the 
indices of the electronic states, and for each electronic state, one sums over the spectrum of eigenstates of the ionic 
system in that electronic state. We have then that there is a contribution to the sum in Eq. 1 from each of the the 
excited electronic states, 5„(g, w) of the form 

(10) 

We are going to approximate the electronic matrix element in Eq. 10 by evaluating it at the equilibrium positions 
of the heavy particles. Since only ionic states for which the struck particle has a momentum in the vicinity of hq 

will have any significant overlap with the state < <I>o.o|e we can assume the \^o.j > that are important in the 

sum have this property. Then this must be true of the states \^n,i > s-s well. Furthermore, for the heavy particles 
which have not been struck, henceforth referred to as the bath, the wavefunctions and energies in both the ground 
and excited states are unaffected, initially, by the scattering event. Hence the change in the heavy particle energies in 

both the excited and ground electronic states is dominated by the kinetic energy of the struck particle, ^5^- Thus 

we can replace {En^i — Eo,j) in the denominator of Eq. 10 by {E^ — Eq). The electronic matrix element is then 
independent of the heavy particle energy indices, and we can use any basis for the particle wavefunctions that we 
like. In particular, we choose a plane wave basis. This also allows us to evaluate the action of Pi , which we will 
replace with Uq. This approximation neglects the momentum of the particles in the initial state compared to hq when 
evaluating the matrix elements. In the energy delta function, we will keep the full kinetic energy ^^^^1 in the 
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final state, and replace _Eo.o by Eq + j^. None of these approximations are essential for our main point, but serve to 
simplify the result. We find then that the excited electronic states contribute to S{q,u)) 

Er^SM^) = S„|^^. < ao| ^^,^^,^, K > P / nim^^ (< - 4)/^ - - ^)dp (11) 

where n{p) is the diagonal matric element of the one particle density matrix for particle 1 in the momentum 
representation. We see that these terms give the same functional form in q and lu as the usual impulse approximation 
result but shifted to high energies by the difference in energies of the electronic states. 

The total intensity in these terms is compensated for in the sum in Eq. 1 by the reduction in intensity of the 
scattering in which there is no change in the electronic wave function. The ground state of the electronic system, in 
the presence of the fast proton, is no longer |ao >, but is given by Eq. 9. When used in the sum in Eq. 1, it must 
be normalized. The normalization of this state leads to a reduction of intensity of the primary signal, again with the 
same approximations, by an amount AI 

This just compensates for the intensity transferred to the excited electronic states. To the extent that one measures 
only the intensity in the primary peak, Eq. 12 gives the apparent deficit in the crossection of the struck particle. We 
sec that this varies with and 1/Af^. as well as the particulars of the electronic system. 

Although a detailed calculation is necessary to evaluate the magnitude of the deficit in the intensity of the main 
peak, we can get some idea of the size of the effect by introducing dimensionless variables. A characteristic length 
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scale for molecular systems is the radius of the Bohr orbit, a. A characteristic energy scale is Since 

^ = (13) 

where the sum in Eq.l3 is over all the electrons in the system, we find that an estimate of the magnitude of the effect 
is given by 

2Neffh^q/a 

where iVe// is the effective number of terms in the sum appearing in Eq. 13. 

For a q of 50A~^ and A^e// of 3, AI « .15. Using a perhaps more realistic scale for water, which would be the 
energy gap of about 6 ev, we would have an increase by a factor of four. The result is also very sensitive to the value of 
Neff. It is clear that without evaluating the integrals, there is considerable uncertainty in the magnitude of the effect, 
but that nevertheless, it is potentially of the right order of magnitude to explain the deficits. However, the result 
above cannot explain two features of the phenomenology of the deficit as it is observed in a variety of systems. The 
deficit is not observed to vary with [3,9] in the Neutron Compton Scattering measurements on water. Although 
the results of the Electron Compton Scattering experiments in formvar [2] are consistent with a dependence of 
the deficit for the small range of q measured, the great bulk of the available data, over a much wider range in q, is 
obtained using neutrons. The data on the deficit in Nb and Pd hydrides [lOjdoes show a strong q dependence, which, 
however saturates at large q. The saturation value of AI is about .4, which is four times the value in the deuterated 
Nb hydride, as the theory above would predict, but the deuterated material shows no variation with q. 

Also, as it stands, there is no explanation for the dependence of the deficit on the relative concentration of deuterium 
and ordinary water. [1] Since the matrix elements and the density of electronic states that determine the amplitude 
of the effect are presumably nearly the same in D2O and II2O, we would expect no concentration dependence from 
the results above. Both of these inadequacies point to the coupling of the electronic states in the presence of a fast 
proton being too strong to use the weak coupling perturbation theory presented so far. The efii'ect of the higher order 
terms in the coupling may be accounted for by introducing linewidths for the excited states due to transfer of energy 
in the electronic system to the ions. These linewidths are strongly enhanced by the presence of the fast proton. 

The time scale for the relaxation of the density correlation function, (^^)~^ [8], where Ap is the momentum 
distribution width, is about 1 femtosecond in water, for typical momentum transfers of 35^1^^. This corresponds to 
a momentum width for the protons of about 5A~^ [11]. It is known that the relaxation time for the lowest excited 
electronic state in water to the so called solvated electron state, is about 150 femtoseconds, [12] considerably longer 
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than the time scale we arc interested in. However, that is when the protons are moving at velocities comparable to 
above. When one of these protons has momentum of 30-40 A^^ the relaxation rate could be a factor of fifty faster, 
since this rate has a component that varies as q^. That rate is comparable to the typical time scale of the scattering, 
and the decay of the excited state can play a major role in determining the intensity deficit. Furthermore, it is known 
[13] that the excited electron wave function extends over 5-10 water molecules, and that the rate of transfer of energy 
to the protons or deuterons depends on a decoherence time scale that is longer by a factor of two in D2O than in H2O, 
providing a mechanism for the concentration dependence of the deficit. Although we are not be able to compute the 
electronic properties needed to make this argument precise, it is plausible that a stong coupling theory that includes 
the decay of the excited states can account for the concentration dependence in water. 

To go beyond lowest order perturbation theory, we will use a mode coupling approximation to include the decay of 
the excited states. The intermediate scattering function, S(q,t), is 

S{q,t) =< 0|e-*«-^''^e*«-^i^(*)|0 >=< 0\e-'^-'^^e^^"°+^^'>*e'^-^'e-*'^"°+^^*\0 > (15) 

Since 

s|^$o,o(^) = (Eofi - KiR))^o,oiR) (16) 

it follows that the action of the Hamiltonian on the ground state in Eq. 15 can be replaced by 

■ ' ■ 

g-tBo.ot^P = e-^^o(^)*e-*^5*'*$o o(ij) (17) 

Having done this, the function e^s"'^o(^)* commutes with the density operator, e^^-^-^ ^ and can be combined with the 
Hamiltonian acting on the final state. Tlicax; is an additional term, arising from the commutator of the kinetic energy 
term in Eq. 15 with Eq(R). This gives the effects of the forces acting on the heavy particles. This term has no direct 
effect on the impulse approximation limit, [8], g ^ 00, t — + 0, qt finite, as it is of order qt^. It leads to traditional 
final state corrections of order 1/q that are included in the data analysis. Thus we can include the ground state 
electronic energy in the exponent in Eq. 15. The kinetic energy terms that do not involve the struck particle can 
also be approximated as commuting with the electronic energy, since this is the essence of the Born-Oppenheimer 
approximation. The noncommutative terms are responsible for the thermalization of the electronic states. These lead 
to the relaxation processes mentioned earlier, and occur on a much longer time scale than the one we are interested 
in. We will therefore neglect the commutator of the kinetic energy of the unstruck particles with the electronic 
Hamiltonian. We can then commute the kinetic energy for the unstruck particles through the density operator and 
eliminate it from the direct time evolution of the density operator. 
Finally, we can approximate the operator Pi^ by 

P^^ = Pl^ ^ 2Pl .Pf (18) 

where the operator Pi acts only on the proton wave function, and the operator Pi acts only on the parameters R 
of the electronic wavefunction. In its action on the ground state, we will ignore the second term on the right in Eq. 
18, which is, again, the Born Oppenheimer approximation. 

We will also replace P^ by hq in the second term of Eq. 18, when this operator acts on final states for which the 
momentum is approximately hq. The full Hamiltonian acting on the final state is now the sum of two commuting 
terms, allowing the operator Pi to act on e'''^^ we can then express the intermediate scattering function as 

S{q, t) =< 0|e-*«-^''ie^(^''-®o+^'*e'«-(^''i+^*)|0 > (19) 

If we neglect the perturbation — p-^ in Eq. 19 we obtain immediately the impulse approximation result for S(q,t), 
making clear again that it is the mixing of the Born-Oppenhcimer cigcnstatcs that is the source of the intensity deficit. 

Eq. 19 can be expressed in more detail, making use of the identity e"2^*e"^'''*'^e*^''^"M )* = e*'J (^i+"M )*, as 

S{q,t) = e^*y" dfd^$S,o('R)ao(^^)e^^^°"-^°^^^+^^*e'^tao(f,^)$o,o(^) (20) 

For each value of R, we see that we now have a standard perturbation problem for the action of the Hamiltonian 
on the electronic wave function. 
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It is clear from the above, that it is the evolution of the electronic state in the presence of the fast proton that must 
be calculated to obtain the deficit. We choose to do the perturbation calculation for the electronic propagator in the 
frequency domain rather than directly in the time domain. Using standard projection operator methods for the the 
resolvent operator for the electronic Hamiltonian we can show 

< ao\[z -Ho + E,{R) + ^I'Vo >= [z - Toiz)]-' (21) 

where 

ToW - KA^f < «o U/^\,J a„ X a„|[z - {Ho - E,)/n - ^t'^m >< a„^|-7^^]->o > (22) 



The prime on the summation in Eq, 22 indicates the that the the ground state is omitted from both sums. Fur- 

thermore, the perturbation operator should be understood as having no matrix elements with the ground 

state. 

We will approximate the sum in Eq. 22 by keeping only the diagonal terms. These are the leading terms in the 
perturbation series in powers of q. In this case, we have 

^o{z) = S„ ^ , , (23) 



where 



rn(^) = S™ (T,' l"^Jn r ( ^ (24) 

The prime on the sum in Eq. 23 indicates that the ground state is omitted, and in Eq. 24 that both the ground state 
and the state m=n are omitted. Eqs. 23,24 are an approximation, which should be understood as self-consistently 
defining the self energies. We will not attempt to calculate these. For our purposes, it is sufficient to note that these 
self energies are proportional to for small q, and can become independent of q for large q. 

The energies in Eq.23 depend on the position of the proton, which is to be averaged over in Eq. 20. The effect of 
this averaging is to produce a distribution of energies in the resolvent expression Eq. 23 that is independent of q and 
will therefore be neglected in comparison with the self energy of the excited states. However, while formally smaller 
in powers of q, the electronic energies are large compared to and their variation can introduce a time scale 

comparable to the decay time of S{q, t) for moderate values of q. In the work of Schwartz et al [13], it is shown that the 

decoherence factor < es'(^'>(^)~^o(^))* > that arises from this averaging has a decay time of about 5 femtoseconds for 
deuterium, and about half that for protons. [13] Therefore the numerical value for ro(^) can be strongly influenced by 
this averaging, which must also include, in a realistic calculation, the average over spatial configurations of surrounding 
molecules. We believe this difference in decoherence time is responsible for the concentration dependence of the deficit 
in H2O-D2O mixtures. The shorter decoherence time produces smaller relaxation rates, smaller values for ro(2;) and 
hence larger deficits. We neglect it here for simplicity only, in order to focus on the main point, the qualitative 
behavior of the intensity deficit with q, and evaluate the energies at the equilibrium position of the proton. The 
electronic propagator(resolvent) then becomes independent of i?i and the momentum of the proton is conserved by 
the dynamics. We find then, that the Laplace Transform of Eq. 19 is 

Ma 



/ dte-"'S{q,t)dt = -i / dpn{p)[z - ^ - ^ - S„ ■ , T" ^ (25) 

Jo J 2M M " z - (E' - E')/h- ^ - ^ -r (z - ^ - ^) 



We can recover the results of the previous weak coupling analysis by ignoring Tn{z) in the denominator of the 
expression in the sum in Eq. 25. The neutron scattering function, S{q, w) in Eq. 1 is obtained by letting z— > w — ie 

in Eq. 25 and taking the real part of the resultant expression. The amplitude of the pole at a; = ^ + ^ is 

^/(^ _ '^'""j^) \z=ui) where ro(z) is the sum in Eq.25. This is the same result as Eq.l2, to lowest order in q. With the 
inclusion of linewidths for the excited states, and hence lineshifts as well, we have the possibility of this intensity deficit 
saturating as q^ 00. What is required is that the lineshifts become comparable to the energy differences from the 
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ground state for the excited states for which the matrix elements in the numerator are significant. While we cannot 
show that this occurs for the values of q in the experiments, since we cannot evaluate the electronic matrix elements, 
we can show that in the limit of very strong coupling, in which the perturbation matrix elements are comparable or 
greater than the energy level separations, such a result holds. In Eq.20, we can write 

e^(^^o--So(«)+T-)t « e*^"°-^oiii))te—^e^'^"^ (26) 

if we neglect higher order terms in the Baker-CampbcU-Hausdorf formula. The last term in Eq. 26 is of order qt^. 
While this is formally of higher order than the terms kept in the impulse approximation limit, the coefficient is 
determined by the forces on the proton in the excited states, and is on the scale of electronic energies. For times of 
the order of the characteristic times that are actually achieved in the experiments, the last phase factor in Eq. 28 is 
approximately 

(qAp\-l 

^ ' « 9 (27) 

Ap ^ ' 

and cannot be neglected. It has been included in the discussion above. It is nevertheless useful to consider the very 
large q limit, where this phase factor can be neglected. We then find that S(q,t) can be written simply as 

S{q,t) « ewt| dmi{R)MR+^) < ao{R)\ao{R + ^) > (28) 

The notation q{ in Eq. 28 indicates that it is a vector corresponding to a displacement of the coordinate of particle 1. 
It should be noted that the electronic matrix element is independent of R at t=0, and will be weakly dependent on R 
for subsequent times, while the nuclear wavefunction depends strongly on i?i Ignoring the position dependence of the 
electronic overlap, we can then write the nuclear overlap integral in the momentum representation, and obtain the 
usual impulse approximation result for S{q, t) multipled by a time dependent factor. Our point here is that this factor 
approaches a constant in the large q limit for times of the order of the characteristic time. Since the time for which 
S{q,t) is significant is (^^)~^, the electronic overlap integral, which is responsible for the intensity deficit, involves 
the two states separated by a finite amount,-^. Or using the uncertainty relation, by a separation comparable to 
the localization distance of the proton in the ground state. The reduction in intensity, therefore, is independent of q 
in the very strong coupling limit, and proportional to q^ in the weak coupling limit. While we are not in the very 
strong coupling limit, we think it plausible that the experiments are being done in an intermediate region, with the 
perturbation comparable to energy level differences, and showing a weak q dependence. This could lead to a variety 
of dependences of the deficit on the electronic properties, as observed. While most experiments show an increase in 
the deficit with q, for instance, recent experiments on LaH2 and LaH3 show a slight decrease with increasing q at 
large q. [14] 

If the intensity is transferred, it is important to know why it has not so far been observed. In water, the minimum 
excitation energy is about 6 ev, so that the intensity that is shifted by exciting the electronic system would be on the 
high energy side of the main signal by at least that amount. It would continue to higher energies by an amount that 
depended upon the sum of the widths of the occupied and unoccupied levels that were involved in the transitions. 
Using recent calculations of Cabral do Couto ct al [15], we find this to be on the order of 30ev if the dominant 
contribution to the sums in Eq.23 come from the valence and conduction band, but the widths could be much larger 
if there is significant excitation of higher and/or deeper bands. It is diflacult to say where the center of gravity of the 
shifted spectral density would lie, but it is easily conceivable that it lies at energies that correspond to transit times 
at Vesuvio of less than 50 milliseconds, which is where the data for water stops. It is also the case, that because 
of the collapse of the phase space for short times, where a large range of momentum transfers are incorporated in a 
small change in the time of fiight, the signal may be missed due to detector saturation. The data for LaH2 and LaHS 
is more problematical, as there is a band gap of about .5 ev in LaHS, while LaH2 is a metal [16] and yet there is 
very little difference in the magnitude of the intensity deficit. This could be explained by the center of gravity of the 
absorbtion being at energies much higher than the gap. 

It is evident that to give a quantitative explanation for the intensity deficit in various materials, it is necessary to 
be able to calculate the integrals over the electronic levels with some accuracy. We have not done this, and indeed, 
it is possible that the effects we are calculating may be too small to explain the observations. What we have shown 
is that a straightforward physical effect, the mixing of Born-Oppenheimer levels and an increase in the linewidth of 
excited electronic levels when a proton of large momentum is present, effects that must be operative at some level in 
any case, have the capacity to account for the phenomenology of the intensity deficits in Neutron Compton Scattering, 
and we have provided a framework in which to accurately calculate these effects. 
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